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In this work, a multichannel Kondo-lattice model is studied in the ther- 
modynamic limit. The conduction band is described by a constant hopping 
amplitude between any pair of lattice sites. For this system we have obtained 
the exact thermodynamical properties and the ground-state energies. In the 
limit of strong interaction between the conduction electrons and the impurity 
spins, the wavefunctions take the Jastrow product form. 





1 Introduction 



The Kondo-lattice model consists of one or several electronic conduction 
bands (channels) interacting with some impurities through a spin exchange 
term [I], |2|]. This model is relevant for the studies of the so-called Kondo- 
insulators and heavy fermion rare-earth compounds. There are considerable 
recent interests in the study of magnetic impurity effects in conduction elec- 
trons [3-7]. 

For one dimensional Kondo-lattice, various numerical and theoretical in- 
vestigations were carried out [8-13]. It was found that at half-filling any arbi- 
trarily small interaction would make the system an insulator [8-13]. One in- 
teresting open question is whether dimensionality plays a role in this Kondo- 
insulator, similar to the Mott-insulator due to electron-electron correlation. 

In this work, we consider a Kondo-lattice model with multichannel con- 
duction bands. In the extreme situation, where the conduction electrons have 
constant hopping amplitude between any pair of lattice sites, we show that 
one can rigorously construct the free energy of the system in the thermo- 



dynamic limit. This generalizes the previous result [l4j to the multichannel 
case. In particular, at zero temperature, our exact solutions demonstrate 
interesting metal-insulator phase transitions in the system when the filling 
numbers of the electrons vary. In the limit of infinite interaction, the ground- 
states take Jastrow product form. It remains unclear how to define effective 
masses of electrons in this simplified model, and how to study heavy mass for 
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the electrons induced by the impurity spins. Nevertheless, in this simple case, 
one can write the exact solutions for the ground-state and the free energy, 
and we report these results and their derivations in the following sections. 

2 The multichannel Kondo-lattice model 

The hamiltonian of the M-channel Kondo model on a lattice with exactly 
one impurity spin at each site and an arbitrary number of orbital conduction 
electrons has the form 

H = H + JV (1) 

where H is the kinetic energy of the electrons, J is the coupling constant, 
and the spin-spin interaction term V is given by 

L M 

^EEl-t (2) 

i=l m=l 

with S{ the spin operator for the impurity at the site i, L the number of 

— * 

sites, m labels the channels, and Sf m the spin operator for the electron at 
the site % in the m-th band. 

Introducing the creation and annihilation fermionic operators (cf am , c iam ) 
for an electron at site i, in the m-th band, with spin a, we have 

Sim = 9 C tam^afiCif3 m (3) 

Z a,/J=U 

where a are the Pauli matrices. Similarly we introduce creation and annihi- 
lation fermionic operators fi a ) for an impurity at site % with spin a and 
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express the spin of the impurity at site i as 

S( = \ E f&tafifit,. (4) 

Clearly the c-operators commute with the /-operators. 

Let H denote the Hilbert space with an arbitrary (0, 1, or 2) number of 
impurities at each site. The condition that there is exactly one impurity at 
each site means that we restrict ourselves to the subspace H where 

E = for all i = l,...,L. (5) 

In our model we consider a kinetic energy if with constant hopping 
amplitude t, i.e. 

M 

H = ~t E EE c £m<W (6) 
l<^j<if=T,i rn=l 

Because of the special form of H , the dimensionality of the lattice is irrel- 
evant and the system is basically one dimensional. In the next section we 
shall thus discuss the one-dimensional system in the thermodynamic limit 
L — > oo. 

Let us introduce the Fourier transform of the electronic operators by 

4am = -jz £tf=i eikjc jv™ and c kam = ^ Ej=i e'^ ''c jam where fc is in the 
first Brillouin zone. In terms of these operators, the kinetic part H of the 
hamiltonian is diagonal 

M 

H o = E EE #) c L c ^ ( 7 ) 

fceFBZcr=t,J, m=l 

with the dispersion relation 

e(k) = -t'LS ok + t (8) 



and t' = t. However, in the following, we will consider t' and t as two 
independent parameters. 



3 Cluster expansion 

Although we are not able to diagonalize the hamiltonian (|I|) for the finite 
system, we can, as in |14| or ||I5| , use a cluster expansion to obtain the 
grand-canonical potential Q/L in the thermodynamic limit. The potential Q 
is defined by 

Q = -^-\nZ (9) 
with Z the grand-canonical partition function 

Z = Tr n e-? K (10) 

where 

M 

K = V™N e m = K + JV 

m=l 
M 

K = H -J2^mN e m . (11) 

m=l 

fi m and are the chemical potential and the number of electrons operator 
of the m-th band. 

In the definition of Z, the trace is performed over the Hilbert space 7i 
where at each site there is exactly one impurity. It is however more convenient 



to work in the Hilbert space 7i. For this, following ||16||, we introduce the 
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impurity isospin at site i by 



r( = \ E f&apfo (12) 

Z a,/3=T4 



where /j is the Nambu spinor: 



fi 



(13) 



fn 
\f& 

Let us consider the sum of the spin and of the isospin at site i: S{ + rf . 
In the subspace where at site % we have exactly one impurity, S{ acts as a 
1/2 spin operator and f( as the zero operator. In the subspace where at 
site i we have zero or two impurities, S{ acts as the zero operator and t( as 
a 1/2 spin operator (the zero impurity spin state corresponding to a down 
spin and the two impurities state corresponding to a up spin). Thus, if we 
replace in K the spin by the sum of the spin and of the isospin, we have 2 L 
different subspaces in which this new operator acts similarly. This allows us 
to rewrite the partition function with a trace on the whole Hilbert space: 

Z = Tr n e-P K ^ = L TrH e-^ §f+7l \ (14) 

Instead of using the expressions (|j) and (|I^) for the definitions of the 
spin and of the isospin, we prefer to use Majorana fermions. For this, we 
decompose the /-spinor in real and imaginary part by: 

/ 




1 / -vl + ir if 



V2 



r)f + irft 



(15) 



where 77°, r/j, rjf and rjf are the Majorana fermions. These operators are 
hermitian and satisfy anticommutation relations: {rjf, rfy} = 5ij5 a b- We define 
their Fourier transform 

% a = ^E^ a (16) 

which satisfy the anticommutation relations {77^, rj^,} = Sk-k'^ab- In terms of 
the Majorana fermions the sum of the spin and of the isospin has a simple 
form: 

Si + ri = - l -ff i Afi i (17) 

where ffi = (r/j, rjf, rjf). Using the commutation relations of the Pauli matrices 
(a — — (i/2)a A a), we can rewrite the spin exchange term as 

Y L M 

V = -oJ2 E E^m^A^feA^). (18) 

i=l a,/3=1,lm=l 

Let us now apply the cluster expansion method to the thermodynamic 
potential: 

1 00 

fi = fio"-:E jnW n (19) 
P n=l 

with Q = -i In {^Tr H e- pK °') and 

W n = (-l) n f dr x r dr 2 ... / T "" 1 dr n <y(r 1 )y(r 2 )-..y(r n ) >g . (20) 
Jo Jo Jo 

V(t) is the free Heisenberg representation of the potential, and < • ■ ■ 
denotes the grand-canonical average of the system without interaction, taken 
only over connected diagrams: 

V(t) = e TKo Ve- rKo (21) 
<V(r 1 )V(r 2 )---V(r n ) > = Tr^poV^V^) ■ ■ ■ V(r n )) (22) 
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where po — e~^ K ° /Trne~ l3K ° is the density matrix of the model without 
interaction. 

Since K is diagonal in the Fourier space, we have for the electronic 
operators 

C kam( T ) = ^ Tc kam (23) 

Cfco-m(T) = 6 ^ ^ ^ 11 ^ Ckam- (24) 

Let us remark that c kam (r) and c kam (r) are no more adjoint to each other. 
Moreover the Majorana fermions commute with Kq and thus they do not 
depend on r: 

Vi(r) = (25) 

This allows us to write explicitly the interaction term V(r) in the Fourier 
space: 

1 _ M 

V(t) = -— HE 4 1 -fc 2 +fe 3 +fc4,oCfc iam (r)c fe2/3m (r)(o ! Aa)^(r/ A; 3A7/ fc4 ). 

Siv k 1 k 2 k 3 k 4 eFBZ a,(3=],lm=l 

(26) 

In order to calculate the grand-canonical averages, we use Wick's theorem, 
which tells us that any average of product of operators can be expressed as a 
sum of product of averages of two operators. We can compute these averages 
explicitly: 

< 4amVk >0 = < C kcn jf k >0= 

< C kam C k'a'm' >0 = < C kam C k ' a' m> >0= 

e -/3(e(fe)- Mm ) 

< C kam C k'a'm' >0 = £fcfc'<W<*mm' J + e _ /3(e(fc) _ Atm) 
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< c kam c k'a'm' >o — Skk'Saa'S mm' 



1 -f e -P{e{k)-Hm) 

< Vtvl' >o = ^5k,-k'5 ab . (27) 

Since e(k) = —t'L5 k + t, these averages do not depend on t' except those 
with momentum zero. However, in the thermodynamic limit, and for if > 0, 
we have 



lim < 4 al (Ti)c 0(7l (Tj) > 

L— >oo 



n = tj 



Ti > Tj 



■ Ti — Tj < (3 

lim < CoriCrOcJrtfo) > = { ' (28) 



L— >oo 



On the other hand, for if < 0, 



lim < c&^rOcarffo) > 

L— >OQ 



Ti ~ Tj = (3. 

T t ~Tj=(3 
Ti - Tj < (3 



1 T' — T ' 

lim < c^y^Tj) > = <! ' (29) 



L— >oo 



Tj > Tj. 



Therefore from Eq. (|20|) , (|28[) , (|29|), the coefficients W n , n > 1, become 
independent of £' in the thermodynamic limit. Assuming we can permute 
the limit L — > oo and the sum over n, we are left with the corresponding 
expression for the model, with t' = 0: 

T / m \ -| oo 

Q = -2MLt'6(t') - - In 2 [] (1 + e^^f --^ J n W n (t' = 0) + 0(1) 

P V m=l / ^ n=l 

1 oo 

= -2MLt'6{t') + Q (*' = 0)--2 J n W n (t' = 0) + 0(1) 

P n=l 

= -2MLt'6{t') + Q(if = 0) + O(l) (30) 



where 6(f) = 1 if f > and 9(f) = if f < 0. = 0) is the thermody- 
namic potential for the system with f = 0. 

We have thus shown that in the thermodynamic limit the potential is 
made of two contributions in which the dependance in f and J separate. 
The first contribution is trivial and we can compute explicitly the second 
one, since it is site by site diagonal. In the next section, this fact will be used 
to obtain the ground-state energy. 



4 Ground-state energy and phase transitions 

In this section, we study the zero temperature properties of the system. The 
ground-state energies will be obtained as functions of electron filling numbers. 
As the filling numbers vary, one finds metal-insulator phase transitions. 

From Eq. fl30|) , the number of electrons in each band and the energy are 
given by 

Wm = _ «L) = _««Lz°>) =NU < = 0) (31) 

M d3tt\ M 

E=Yj VrnN m + = ~2fML6(f) + t £ N m + E m \ (32) 

m=l t' / m um m=l 



where 



L M 



E m *(/5,/x 1 ,...,M = ^EE <SL-Sl > % (33) 

i=l m=l 

and 



< A>i = TrA Pi 
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Pi ~ Tre-f* H * 

M 

Hi = £ [ J SL -Sl + (t- ^m)NL] ■ (34) 

m=l 

We are thus left with a system of spins with only on-site interaction. 

At this point several definitions of the ground-state can be introduced. 
In the grand-canonical ensemble, taking j3 — > oo with . . . , /i M ) fixed, we 
obtain only some specific values for the electron density n. For example for 
a single band and J > 0, we have 

n = for pi <t — ^ 

n = | for fj, — t—^- 

n=l for ^t-Mt+Mf (35) 



n=| for yU = t+ M 

3J 



3 ^ ^ * i 4 

n = 2 for pL> t+^f. 



4 

In the canonical ensemble we can consider the ground-states either for 
fixed (JVi, . . . , iVkf), or for fixed N = X)m=i with respect to the hamilto- 
nian 

M 



if = -2t'LM9(t') +tJ2N m + E i 

m=l 

L M 

E int = JY^Y,SL-S{. (36) 

i=l m=l 

Therefore we only need to compute the ground-state energy of the spin- 
exchange interaction, which we shall now obtain explicitly for one and two 
bands. 
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For a single band, we have to diagonalize the on-site exchange operator 
S% ■ S- . For a given site i, the size of the Hilbert space is 8. The four basis 
states with zero or two electrons on this site have energy zero. For the four 
states with one electron we have a triplet with energy J/4 and a singlet with 
energy —3 J/4. Using these results, we recover immediately the ground-state 
energy obtained in ||14|| , both for the antiferromagnetic (J > 0) and the 
ferromagnetic (J < 0) case, i.e. for the antiferromagnetic case: 

\ -^-N N < L 

Eq S ={ ~ ( 37 ) 

[ ~^{2L-N) N > L, 

for the ferromagnetic case: 

ESs=\ * ~ (38) 

[ {{2L-N) N > L. 

For the two bands case, we have to diagonalize the on-site operator 5!f (1) • 
S{ + Sf(2) ■ S{ . The size of the Hilbert space for the site % is 32. The 8 states 
with either zero or two electrons in each band at site i have energy zero since 
the electronic spins are zero. For the 16 states with one band occupied by 
zero or two electrons and the other band occupied by one electron, we have 
again a triplet with energy J/4 and a singlet with energy —3 J/4, since the 
energy for the band with zero or two electrons is zero. Finally, for the 8 
states with exactly one electron in each band, we have a quadruplet with 
energy J/2, a doublet with energy — J and a doublet with energy zero. The 
ground-state energy Egg for the system with N\ electrons in band 1 and A^2 
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electrons in band 2 can thus be computed explicitly and is given by E a for 
(Ni, N 2 ) in the domain V a , a — 1, ... ,8, represented in figure 1 below. 




(a) (b) 



Figure 1: Ground-state energy for J > (a) and for J < (b). Solid lines 
indicate filling numbers where phase transitions occur. 

For the antiferromagnetic case (J > 0), the ground-state energies are 
given by 

3 



E 1 = -| J(N, + N 2 ) E 5 = -f J(2L — Nx + N 2 ) 

E 2 = -\J{2L + JVi + N 2 ) E & = -\J{AL - N x + N 2 ) 

E 3 = -\ J(6L — Ni — N 2 ) E 7 = -\J{AL + N X - N 2 ) 

E 4 = -f J(4L -N x - N 2 ) E s = -\J{2L + N t - N 2 ) 



(39) 
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and for the ferromagnetic case (J < 0) 

E X = E 2 = \j{N x + N 2 ) E 5 = E 6 = \J(2L -N x + N 2 ) 

E 3 = E 4 = f J(4L -Nt- N 2 ) E 7 = E S = \j{2L + N X - N 2 ). 

(40) 

The differences between the two cases J > and J < can be seen as 
follows. In the domain V>\ where N\ + N 2 < L, the lowest energy states are 
obtained with zero or one electron at each site in any of the two bands, which 
gives Ei = — §JiV if J > and E 1 = |iV if J < 0. On the other hand in 
the domain T> 2 we have (2L — N) sites with one electron and (N — L) sites 
which are doubly occupied with one electron in each band and total electron 
spin 1. In this domain T> 2 , if J > the total electrons plus impurity spin is 
1/2, the three spins form a doublet, and the energy is E 2 = — j(2L + N); in 
this case dE^/dN is discontinuous at N = L. However if J < 0, the total 
electrons plus impurity spin is 3/2, the three spins form a quadruplet, and 
the energy is E 2 = \ JN] in this case there is no discontinuity of dE^g/dN 
at iV = L. It is interesting to note that even though the two electronic spins 
are not directly coupled in the hamiltonian, they do interact indirectly via 
the impurity spin. For a system with Ni + N 2 = N fixed, it is easily seen that 
the ground-state is obtained for (N\, N 2 ) in the domains V a , a = 1, 2, 3, 4. 

The exact expressions for the ground-state energies allow us to study the 
conducting properties of the system at zero temperature. For the single band 
case, following an idea due to Mattis [T7|, 0], we define /i + = E GS (N + 1) — 



Egs(N), /I" = E GS (N) - E GS (N - 1) and A/x = //+ - \T . If A/x = 0, 
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the system is conducting; if A/i > 0, the system is insulating due to the 
apparition of a gap A/i. For a single band, it follows from (|3?D and fl38|) that 
A/i = except for iV = L where A/i = 3J/2 if J > and A/i = -J/2 if 
J < 0; therefore, as shown in |I4| , the system is always conducting except 
at the point N = L where it is insulating. For the two bands case, if we 
consider the ground-states with N\ and N 2 fixed, we have to introduce two 
parameters A/ii and A/i 2 referring to the two bands, 

A/i! = EcsiN, + 1, N 2 ) - 2E GS (N h N 2 ) + Eqs^ - 1, N 2 ) 
A/x 2 = E GS (N U N 2 + l)- 2E GS (N U N 2 ) + E GS (N U N 2 - 1). (41) 

The system is insulating if and only if a gap appears in the two bands, i.e. 
if A/i! > and A/i 2 > 0. 

For the antiferromagnetic case, the system is insulating at the point N\ = 
N 2 = L where A/ii = A/i 2 = \ and on the four lines Ni + N 2 = L, Ni — N 2 = 
L, N 2 — Ni = L, Ni + N 2 = 3L where A/i X = A/i 2 = ^. We can thus conclude 
that a metal-insulator transition occurs on the four heavy lines of figure 
1(a) and at the point Ni = N 2 = L, as we turn on the antiferromagnetic 
interaction between electrons and impurities. For the ferromagnetic case, the 
system is insulating only at the point Ni = N 2 = L where A/i! = = ~'\ 
and thus at this point a metal-insulator transition occurs as the ferromagnetic 
interaction is introduced. Similarly, if we consider the ground-states with 
Ni + N 2 = N fixed, then for J > the system is insulating if iV = L, 2L, 3L, 
while for J < it is insulating only for iV = 2L. 
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5 Ground-state wavefunctions at J = oo 



In this section, we consider a finite system A with L sites and an arbitrary 
number M of electronic bands. We want to obtain the ground-state wave- 
functions for the case t < and in the limit J = oo. 

If the total number of electrons is smaller than the number of sites, i.e. 

M 

N=Y,N m <L, (42) 

m=l 

then each electron will attempt to form a singlet of energy —3 J/4 with an 
impurity. Therefore in the limit J = oo we can reduce the Hilbert space 
to the subspace where at each site there is either a singlet or an unpaired 
impurity spin. A general basis for this subspace is given by the following 

vectors 



I1/^J0> (43) 



I a > = | Xi, . . . , Xm', & > 

M \ I ■-- x 

= IT IT ~?K^ C x'\mfx[ ~ C x[mfx y \) 

m=l \_x£X m J n =l 

where the N m singlets with electron in the band m are located at the positions 
X m = (xi m , . . . , XN m m), X m fl X m i = 0, and the unpaired impurity spins a = 
(<7i, . . . , <7l-n) are located at the positions Y — A \ U m X m = (yi, . . . , i/l-n) 
with y 1 <y 2 < ... < y L _ N . 

With P the projector onto this J = oo subspace, the projected hamilto- 
nian has the form 

3 7 

PHP = PH P --N (44) 
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where Hq is the kinetic term. Moreover, we can then forget the (infinite) con- 
stant and consider only the term PHqP. To find the ground-state wavefunc- 
tions, we identify the singlets as spinless bosons and the unpaired impurities 
as spin 1/2 fermions. Let us then introduce the hamiltonian 

L M 



t 



h = -o E EE P G F ia Bt m F+B ]m P ( 



G 



(45) 



v£j=l cr=T,i m=l 

where the B fields are bosonic, and the F fields are fermionic. The B fields 
commute with the F fields and J2i=i Bt m B irn = N m . The Gutzwiller projector 
Pg projects on the subspace where at each site there is exactly one particle, 
either one fermion or one boson, i.e. 



M 



-BimBim + F i(J Fi a 1. 



(46) 



m=l 



For this system, the basis vectors can be taken as follows 



a >- 



M 

n 

m=l 



n b; 



L-N 



UK, 1 > 



(47) 



n=l 



The two systems defined respectively by the hamiltonians PHqP and h are 
isomorphic, since there is a one to one correspondence | a >— >\ a > between 
the basis vectors and the matrix elements of the two hamiltonians are the 



same 



< | PH P I a >=< p\h\a> . 



(48) 



As for the Hubbard model 0, pQ] , we try to write the hamiltonian as 
a positive definite form to determine the ground-state energy. Using the 
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relations 

PGF^B irn F jr7 Bf m P G = -F jr7 B+ m P G F+B im if % ^ j 

P G Bt m B im P G = F la B+ m P G F+B m (49) 

we can rewrite the hamiltonian h in the following way: 

y. L M M 

h = -o E E E F ia B? m P G F+B jm + tP Q £ iV m P G . (50) 

j,j=l ct=T,| m=l m=l 

The first part of h is positive definite since t < 0, and the second part is a 
constant. 

If TV = Y^m N m < L — 2, the ground-state wavefunctions are labelled by 
the positions Y — (yi, . . . , uq) and the spins o = (a±, . . . , <tq) of L—N—2 = Q 
fermions. These wavefunctions are given by 

M 

i<wy»>=e e «n 

zi,Z2 X!,...,x M m=l 

U m X m U{z 1 ,z 2 }=^\Y 

(51) 

where in Eq. ([)]]) we sum over the positions -21,-22 of the two last fermions 
with opposite spin, and over the positions Xi, . . . , Xm of the bosons in each 
band: 

X m fl X m i = if m 7^ rri 
X m n{z 1 ,z 2 } = 

To prove that ([|l]) is indeed a ground-state, one shows that it is annihilated 
by the first part of h and therefore ( [5"TD is a ground-state wavefunction with 
energy tN. 
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n bi 



n 



n 

n=l 



^ <x I > 



If TV = J2 m Nm = L — 1, the ground-state wavefunctions have only one 
delocalized fermion and are parametrized by the spin a of this fermion. They 
are given by 



M 



n s >= E e n n 

z x 1 ,...,x M m =\ 

Um-f m U{z}=A 



n b; 



o > 



(52) 



with 



X m fl X m i = if m 7^ m! 



x m n {^} 



and the energy is |(L — 1). 



As for the single band case [nj], if we impose periodic boundary con- 
ditions, we can write the ground-state wavefunctions in a Jastrow product 
form. Let us consider a state with N + fermions with up spin localized at Y + = 
(jji, . . . , ), N_ fermions with down spin localized at Y~ = (jrf , . . . , y^_) 
and Q = L — N + — holes localized at X = (xi, . . . , £q). The ground-state 
wavefunctions are given by 



Ai 



^ GS >= E n 



m=l 
B 



L 



■n^;vvtn^ + Tio> 



(53) 



71=1 



with amplitude 



(54) 



where the function d(n) = sm(mr/L) and the quantum numbers and 
m s are integers or half-integers which make sure of the periodic boundary 
conditions and satisfy the following inequalities 

-(AL + Q + 1) <m h <L- -(AL + Q + 1) (55) 

i(iV + + Q + l)< mh - ms + L< L - I(iv + + Q + 1). (56) 
To see that such states are ground-states, one has to verify that they are 



eigenstates of the hamiltonian (|45|) with eigenvalues tN. To establish this 
result, one has first to apply the up-spin part of this hamiltonian and see that 
| ^cs > is an eigenvector with eigenvalue ~iV. To apply the down-spin part 
of the hamiltonian, it is more convenient to write the ground-state in terms 
of holes and up-spins [3T|. In these terms, the amplitude *&(X, Y + ) has the 
same form as ^(X, Y~) except that rrih is replaced by — m s + L/2 and 
m s by L — m s . In this manner, it gives an energy |iV and the total energy 
is tN. 



6 Summary 

The ground-state energy of our Kondo-lattice model was obtained explicitly 
in the thermodynamic limit for one and two electronic bands. From this so- 
lution the insulating or conducting properties of the system were established 
as a function of the number of electrons in each bands. For the ferromagnetic 
system with two electronic bands a metal-insulator phase transition appears 



19 



at Ni = N 2 = L, as the interaction between the electrons and the impurities 
is switched on at zero temperature. For the antiferromagnetic system, new 
metal-insulator phase transitions also appear for other values of Ni and N 2 . 

For the finite size system, we have shown that the ground-state wave- 
functions can be written in the well-known Jastrow product form. It remains 
unclear how to define the effective masses for electrons hopping with an un- 
constrained hopping amplitude and whether the impurity spins induce heavy 
masses for the conduction electrons. 

This work was supported in part by the Swiss National Science Founda- 
tion. 
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